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Abstract 

We come back to the Cauchy integral equations 
occurring in radiative transfer problems posed in finite, 
plane-parallel media with light scattering taken as mono- 
chromatic and isotropic. Their solution is calculated fol- 
lowing the classical scheme where a Cauchy integral 
equation is reduced to a couple of Fredholm integral 
equations. It is expressed in terms of two auxiliary func- 
tions and C_ we introduce in this paper. These func- 
tions show remarkable analytical properties in the com- 
plex plane. They satisfy a simple algebraic relation which 
generalizes the factorization relation of semi-infinite me- 
dia. They are regular in the domain of the Fredholm in- 
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tegral equations they satisfy, and thus can be computed 
accurately. As an illustration, the X- and Y-functions 
are calculated in the whole complex plane, together with 
the extension in this plane of the so-called Sobouti's 
functions 

Key Words: Radiative transfer equation; Finite, 
plane-parallel medium; Isotropic scattering; Cauchy 
integral equations 



I. INTRODUCTION 

In Refs. [1] and [2] (hereafter I and II), we revisited the Cauchy 
integral equations arising in the simplest radiative transfer problems 
posed in plane-parallel geometry, that is in homogeneous and station- 
ary media with monochromatic and isotropic light scattering. The 
literature on this subject essentially dates back to the sixties. I 3-8 ! 
In the present paper, we deal with singular integral equations of the 
form 

a f 1 dv 

T(a,z)X (a,b,z) - -z X (a,b,v) 

2 Jo v- z 

exp(-b/z)^-z f X (a,b, -v) exp(-b/v)-^— = c (z), (1) 
2 Jo v + z 

which are encountered when solving transfer problems in finite slabs' 9 ' 
The parameters a and b are respectively the volumic albedo of the 
medium and its optical thickness at a fixed frequency (0 < a < 1 
and < b < +oo). The dispersion function T(a, z) is defined in C \ 
{-!,+!} by 

We specify that the integrals in Eqs. (1) and (2) are Cauchy prin- 
cipal values when the variable z is in the range ] — !,+![, with no 
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change in notation: as in I and II, the same symbol is used for a 
sectionally analytic function outside and on its cut. 

The co-function on the right-hand side of Eq. (1) is defined and an- 
alytic in C* = C \ {0}, including infinity. As a matter of fact, co(l/ ' z) 
is the Laplace transform, over the finite interval [0,6], of the source 
function So for the direct field within the slab, vizJ 9 ! 



As a result, this source term is bounded on the right at z = 0, which 
means that the limit of cq{z) exists and is finite when z tends to 
in the (stict) right complex half-plane. This remains true for the 
function z — ► co(— z) exp(— b/z), as can be seen by replacing So(r) 
by So (6 — r) in the integrand of Eq. (3). Whence 

lim cq(z) < +oo , lim [cnf— z) exp(— b/z)} < +oo, (4) 
z^0,!R(z)>0 2^0,SR(z)>0 

where denotes the real part of the complex number z. We note 
that the possibility of cq defined and analytic at z = is not ruled 
out. It happens when Sq(t) = S(t), the Dirac distribution at the 
origin. We have then cq{z) = 1 and co(— z) exp(— b/z) = exp(— b/z). 
These expressions are exchanged when So(r) = 5(b — r), i.e., cq(z) = 
exp(—b/z) and co(— z) exp(— b/z) = 1. If So is a function or the Dirac 
distribution at r g]0, b[, cq{z) and co(— z) exp(— b/z) are not defined 
at z = and the limits in (4) are zero. 

The unknown function Xq depends on the given source term cq and 
the parameters a and b. This dependence is made explicit since the 
0-subscript is intended to remind of the link with Co . This subscript 
is removed when dealing with the particular free term cq = 1, the 
associated A-function being the well-known function of finite plane- 
parallel media. The function Ao can be assumed analytic everywhere 
in C*, since Ao(a, b, 1/z) is the finite Laplace transform of the source 
function S for the diffuse field within the slab. I 9 ] 

The second integral term in Eq. (1) vanishes in a half-space (b = 
+00) when $l(z) > 0. This equation reduces to the simpler form 




(3) 
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where c(z) = cq(z). Thus it is of Cauchy type on [0, 1]. We studied 
this kind of equation in I, and gave in II its unique solution analytic 
in the right complex half-plane. From now and subsequently, we refer 
to the equations of I (resp. II) by their original number preceded by 
the roman numeral I (resp. II). The solution (II, 50) to Eq. (5) was 
written for a free term c verifying the three following properties: 
(i) it is defined in C \ { — 1}, except possibly at z = 0, 
(m) it is analytic in the half-plane K(z) > (including infinity) except 
possibly at z = 0, 

(Hi) it is bounded on the right at z = 0, i.e., the limit of c(z) when 
z — > with K(z) > exists and is finite. 

We have expressed this solution in terms of the particular solution 
H = H(a,z), which corresponds to the free term c(z) = 1. This is 
just the well-known H- function for semi-infinite media. 

In a finite space (b < +oo), the two integral terms in Eq. (1) are 
kept, so that if this equation is written in the form (5), the free 
term c(z) depends on the unknown function Xq. This substantial 
complication prevents us from solving Eq. (1) explicitly, as in the 
semi-infinite case. 

The remainder of this article may be summarized as follows: in 
Sec. II, we use the classical transformation from Eq. (1), which is of 
Cauchy type over the line segment [—1, +1], to two Cauchy integral 
equations over [0, 1]. Their solution is expressed in terms of two aux- 
iliary functions depending on the free term cq. In Sec. Ill, we mainly 
deal with the study of the auxiliary functions associated to the par- 
ticular free term cq = 1, which we denote C+ and It is shown that 
these functions have remarkable analytical properties in the complex 
plane. In particular they satisfy a simple algebraic relation which 
generalizes the factorization relation for semi-infinite media. These 
functions allow the calculation of the X- and Y-functions associated 
to finite, plane-parallel media, together with the extension to the 
complex plane of the so-called Sobouti's functions. We come back to 
the general case in Sec. IV, where it is shown that the general solu- 
tion to Eq. (1) can be expressed by means of the functions and 
C- we studied in Sec. III. The latter are thus the basic functions for 
the solution of Eq. (1) with any free term cq. Finally, an algorithm 
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for the computation of ( + and £_ is proposed in Sec. V. It is based 
on the resolution of two uncoupled Fredholm integral equations with 
regular kernels. 



II. PRINCIPLE OF THE SOLUTION OF 
EQUATION (1) 

Due to its second integral term, Eq. (1) is of Cauchy type over 
[— 1,+1] while it was so over [0, 1] in the semi-infinite case. We are 
thus confronted with a first difficulty: the solution to Eq. (1) can't 
be Holder continuous in [—1, +1] since it diverges at zero. This very 
fact prevents us from applying directly the theory of Cauchy inte- 
gral equations. This difficulty is overcome by replacing the solution 
to Eq. (1) over [— 1, +1] by the one over [0,1] to a system of two 
singular integral equations satisfied by the functions Xo(a,b,z) and 
Yo(a, b, z) = Xo(a, b, —z) exp(— b/z). This classical approach was fol- 
lowed by Busbridget 3 ! and Mullikin et al.I 4 ~ 7 ] among others. 

The functions Xq and Yq are related for any z € C* by 

X (a, b, z) = Y (a, b, -z) exp(-b/z), (6) 
Y (a, b, z) = X (a, b, -z) exp(-b/z), (7) 

and we have for z G C \ {±1, 0} 

dv 

z 



dv , . 

co(z), (8) 



V + z 



a f 1 dv 

T(a,z)X (a,b,z) - -z X (a,b,v) 

z Jo v ~ 

a f 1 

+ exp(-b/z)-z / Y (a, b,v) 
1 Jo 

a f 1 dv 
T(a,z)Y (a,b,z) - -z Y (a,b,v) 

Z Jq V — Z 

a f 1 dv 
+ exp(-b/z)-z X (a,b,v) — ■ — = co(-z) exp(-b/z). (9) 

Z Jq V + Z 

The second equation is deduced from the first one by changing z into 
—z in it, then multiplying the resulting equation by exp(— b/z). 
The coupling between Eqs. (8) and (9) can be removed by sub- 
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stracting and adding them. We obtain 

T(a,z){X TY ){a,b, z) 

a f 1 dv 
--z (X T Y )(a,b,v) — — = c ,±(a, b, z), (10) 
Z v — z 

where, for any z£C \ { — 1, 0} 
co,±(a,b,z) = 

a f 1 dv 
co(z) T [cq(-z) {X TY )(a,b,v)— —]exp{-b/z). (11) 
Z Jq v + z 

This expression shows that the limits 

co ±(a, b, +0) = lim c ,±(a,6, z) (12) 

z^0,!R(z)>0 

are given by 

c ,±(a,6,+0) = lim [c (z) =F c (-z) exp(-6/z)] < +cx), (13) 

z^0,5R(z)>0 

the inequality being due to (4). 

Equation (10) is a Cauchy integral equation of the form (5) for 
the functions Xq^Yq. Since its free term Co,±(a, b, z) satisfies the 
properties (i)-(iii) as given in the introduction, its unique solution 
analytic in the right complex half-plane is given by Eq. (II, 50), viz. 

(X T*o)(a,M) = 
H(a, z) | ^co,±(a, b, +0) + Y[-K(z)]fl"(o, -z)c 0: ±(a, b, z) 

z r +icc dz' *\ 

+^-J ioo HiaA/z'MaA-l/z')-^^. (14) 

Here Y is the unit step (or Heaviside) function extended by 1/2 at 0: 
Y(x) is for x < 0, 1/2 for x = and +1 for x > 0. The integral is 
a principal value at infinity. This relation is valid on the imaginary 
axis (denoted hereafter iM), provided that the value Y(0) = 1/2 is 
adopted and the integral is taken in the sense of the Cauchy principal 
value. 

We introduce the functions 

uo,±{a,b,z) = (X t Y )(a,b,z), (15) 
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a f 1 dv 
vo,±(a,b,z) = co(-z) - -z / (X ^Y )(a,b,v) — ; — . (16) 

1 Jq v 



+ z 



The functions uq,± are denned on C* where they are analytic. 
They diverge on the right-hand side of 0. The functions vo± are 
defined in C \ {— 1, 0}. They are sectionally analytic over C \ [—1,0], 
the integral being a Cauchy principal value for z G] — 1,0[. They 
diverge at -1 due to the integral term, and possibly on the right side 
of due to the co(— z)-term. The divergence at -1 results from the 
behavior of the functions (Xq ^Yo)(a, b, v) on the left side of +1: 
they don't vanish there (see Ref. [10], p. 42). 

It follows from Eqs. (6)-(7), (15)-(16) and (10)-(11) that the func- 
tions uo,± and «o,± are related by 

u ,±(a, b, z) = Tuo,±(a, b, -z) exp(-b/z) , (17) 

v 0i ±(a,b,z) = co(-z) - ^z f uo t± (a,b,v)-^—, (18) 

l Jo V + z 

T(a,z)u 0j ±(a,b,z) = v 0> ±(a,b,-z) T v 0> ±(a, b, z) exp(-b/z), (19) 
and Eqs. (14) and (11) become 



Uo,±(a,b,z) = 

H{a, z) | ^co,±(a, b, +0) + Y[-K(z)]H(a, -z)c 0i ±(a, b, z) 



+ i/r^ a ' i/z/K±(a ' 6 '- i/z/) i^}' 



(20) 



c ,±(a, b, z) = cq(z) T vo,±(a, b, z) exp(-b/z). (21) 



Now multiply both sides of Eq. (20) by T(a, z) in order to eliminate 
uo t ±(a, b, z) with the help of Eq. (19), and substitute the expression 
(21) of cq,± in the right-hand side of Eq. (20). Using the factorization 



8 



RUTILY, CHEVALLIER AND BERGEAT 



relation (II, 27) and changing z into —z, one obtains 

H(a, z)v 0: ±(a, b, z) = -co,±(a, 6, +0) + Y[»(z)].ff(a, z)c (-z) 
±Y[-$t(z)]H(a, z)v ,±{a, b, -z) exp(6/z) 

-^L (22) 

dz' 

x[co(l/z')Tv ,±(a,b,l/z') exp{-bz')] ■ 
We introduce the functions Co,± = Co,±( a > b, z) as defined in C\zR 

by 

Co,±(a, &, z) = H(a, z)v 0j ±(a, b, z) 

T Y[-K(s)]fT(a, z)«0,±(a, 6, -z) exp{b/z). (23) 

It follows from Eq. (22) that they satisfy the integral equations 
Co,±(a, 6, z) = -Oo,±(a, 6, +0) + Y[St(z)]H(a, z)cq(-z) 

-^L H{a -- il - ] (24) 

since fo,±(a, I/- 2 ') can be replaced by (£o,±) + (g, 6, l/z')/H(a, 1/z') 
in the integral of Eq. (22), due to Eq. (23). The functions (Co,±) + are 
the limit of the functions £o,± on the right-hand side of the imaginary 
axis. Relation (24) shows that the latter functions are sectionally 
analytic in the complex plane cut along the imaginary axis, with 
limits (£ 0) ±) + and (Co,±) _ on the right- and left-hand sides of this axis 
respectively. We could extend the definition of Co,± to the imaginary 
axis by calculating the integral in Eq. (24) in the sense of the Cauchy 
principal value; however, such an extension does not seem useful in 
the context of this article. We note in addition that the integral in 
Eq. (24) is a principal value at infinity. It can be transformed into 
an integral over [0, 1] using the residue theorem, which implies the 
existence and uniqueness of the functions £o,± as defined by Eq. (24). 
Here we dispense with the details of this demonstration. 
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Replacing vo t ±(a,b,—z) by Co,±(a> b, — z)/H(a, — z) [9ft (z) < 0] on 
the right-hand side of Eq. (23), we derive the following solution to 
Eq. (22) in C\{iMU{-l}}: 

, , ± (a, 5, z) = g(a>g) ± Y[-K(s)] g(a _ g) exp(6/z). (25) 

Hence ito,± in C \ iR using Eq. (19) and the factorization relation 
(II, 27) 

u 0i ±(a, b, z) = Y[9ft(z)]i7(a, z)( 0i ±(a, b, -z) 

TY[-»(zp(a, -*)Co,±(a, 6, z) exp(-tyz). (26) 

These expressions can be extended to the imaginary axis by letting 
Dft(z) — ► 0* in them and taking into account the continuity of the 
functions uo,± and vo t ± on the iR-axis. We derive their restriction on 
this axis in terms of the limits of the functions £o,± on both sides of 
the axis. Hence for any zq £ iM* 

/ . x (Co,±) + (a,Mo) , 07A 
^(a,Mo) = g(fl>zb) , (27) 

= (Co,,)-(a,Mo) ± (to,±na,b,-zo) 
H(a,z ) H(a,-z ) 



w ,±(a,Mo) = H(a,z )((o,±) (a,6,-z ), (29) 
= TH(a, -z )(( 0) ±y (a,b, z ) exp(-b/ z ). (30) 

We may " invert" these relations by writing on C \ iM. 

Co,±(a, b, z) = Y[3ft(z)]ii"(a, z)v ,±(a, b, z) 

which shows that the functions Co,± give information about i>o,± from 
their values in the right complex half-plane, and about uo,± from 
their values in the left one. 

The solution Xq to Eq. (1) and related functions appearing in the 
left-hand side of this equation are easily deduced from the functions 
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ito,± an d vq : ±, since we have from Eqs. (15)-(16) 




o 




Finally, the only thing we have to do is to solve the problem 
(24), since its solution £o,± leads to the functions uo,± and vo t ± via 
Eqs. (25)-(30), then to the solution to Eq. (1) via Eqs. (32)-(35). In 
the next section, we first solve Eq. (24) for the particular free term 
co(z) = 1. 

III. THE c = 1 CASE: CALCULATION OF THE 
X- AND Y-FUNCTIONS 

We treat the case cq = 1 in this section. The associated functions 
will be denoted by the same symbol as in the general case, just 
removing the 0-subscript. Equation (24) simplifies for co = 1, since 
we have c±(a, b, +0) = 1 from Eq. (12) and 



owing to the residue theorem applied to an obvious contour in the 
left complex half-plane. The functions C± are thus defined on C\iM* 
by the integral equations 




(36) 




xexp(-bz')((±) + (a,bA/z') 



dz' 



(37) 



1+zz' 
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These equations have a sense at z = and yield (±(a,b,0) = 1, 
since the identity (36) is valid at z = [Y(0) = 1/2]. In Sec. V, it 
will be seen that Eqs. (37) can be transformed into two uncoupled 
Fredholm integral equations over [0,1]. A detailed study of these 
equations will be published in a forthcoming paper. It can be shown 
that they admit a unique solution in their domain, which implies the 
existence and uniqueness of the functions (± as defined on C \ iR* 
by Eqs. (37). 

The functions (± are sectionally analytic in the complex plane cut 
along the imaginary axis, with (C±) + and (C±)~ as limits on both 
sides of this axis. These limits satisfy the relations (27)-(30) with 
Co = 1. Comparing the right-hand sides of these equations (with — zq 
instead of z$), we obtain for any z$ 6 iM* 

H(a,z )((±) + (a,b-z ) = H(a,z )((±y(a,b-z ) 

±H(a,-z )(C±) + (a,b,z ) exp(-6/z ), (38) 

H{a,Zo)(C±)~(a,b-Zo) = TH(a-Zo)(C±)~(a,b,Zo) exp(-b/z ). (39) 

Note that Eq. (38) could have been deduced from Eq. (37) and 
the Plemelj's formulae on Cauchy-type integrals J 10 ' 

The functions (± have remarkable analytical properties, including 
the following one, valid over C\iR* 

-[C+(a, b, z)C-(a, b, -z) + O, b, z)(+(a, b, -z)} = 1, (40) 

which solely derives from the definition (37). The proof is given in the 
Appendix. This relation generalizes, in a finite slab, the factorization 
relation specific to semi- infinite media (see the conclusion). It leads 
to the following three properties of the functions u±, v± as defined 
by Eqs. (15)-(16) without the 0-subscripts and with cq = 1: 

-[u-(a, b, z)v+(a, b, z) + u+(a, b, z)v-(a, b, z)] = 1 (z ^ -1, 0), (41) 

— [u-(a, b, z)v+(a, b, —z) — u+(a, b, z)v-(a, b, —z)\ = 

exp(-b/z) (z^0,+l), (42) 
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-[v+(a,b,z)v-(a,b, -z) + v-(a, b, z)v + {a, b, -z)] 



T(a,z) (z^±l). (43) 



To summarize, the Co = 1 case is characterized by the following 
four points: 

1) The X- and Y-functions are the unique solution analytic in C* to 
the singular integral equations (8)- (9) with co = 1. They are thus the 
classical functions of finite plane-parallel media, first introduced by 
AmbartsumianI 11 ! and readily studied by ChandrasekharJ 12 ! These 
functions are connected by the relations (6)-(7) without the 0-sub- 
scripts, viz. (z £ C*) 



T(a,z)X(a,b,z) =1- U(a,b,-z)- fr{a,b,z)exp(-b/z), (48) 
T(a,z)Y(a,b,z) = [l-fr(a, b, z)]exp(-b/z)- fr(a, b ,-z), (49) 

for any z in C\{0,±1}. 



2) The functions X, Y, and can be expressed in terms of 
the functions u± and v± by means of Eqs. (32)-(35), removing the 
subscripts and putting Co = 1 in Eq. (34), which is then valid at 
z = 0. Hence 



X(a, b, z) = Y(a, b, —z) exp(— b/z), 
Y(a, b, z) = X(a, b, —z) exp(— b/z). 

Introducing the so-called Sobouti's functions' 13 ! 



(44) 
(45) 




(46) 



(47) 



the X- and F-singular equations (8)-(9) read 




(50) 



(51) 




(52) 
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^zf 1 Y{a l b,v)-^- = hv + -v-){a l b,z) (z€C\{-l». (53) 
Z Jq v-\-z z 

3) The functions u± and v± follow from (± using Eqs. (25)-(26) 
on C\iR and Eqs. (27)-(30) on iR* . These functions satisfy the 
general equations (17)-(19), which are equivalent to Eqs. (44)-(49). 
When Co = 1, they also verify the properties (41)-(43) resulting from 
the relation (40) between ( + and These three equations may 
be rewritten in terms of the functions X, Y, £x and £y using Eqs. 
(50)-(53), viz. 

X(a, b, z)[l-i x {a, b, z)] + Y(a, b, z^ Y (a, b, z) = 1 (z ^ 0, -1), 

(54) 

X(a,b,z)£ Y (a,b,-z) + Y(a,b,z)[l - £ x {a,b,-z)} = 

exp(-b/z) (z^0,+l), (55) 

[1 - £x(a,b,z)][l - £ x (a,b,-z)] - fr(a,b,z)£ Y (a,b,-z) = 

T(a,z) (z^±l). (56) 

These equations are well-known: (54)-(55) are just the regular X- 
and F-equationsl 11 ' 12 ! and (56) was stated half a century ago. I 14 ! We 
have retrieved the main functional properties of the functions X, 
Y, £x and £y solely from the singular X- and Y-equations (48)- 
(49). This approach is opposite to the classical one, with central role 
played by the regular X- and Y-equations (54)-(55).[ 12,15 J 

4) Finally, the functions (± are calculated from their definition (37) , 
transforming the integral over the imaginary axis by the method of 
residues. Two Fredholm integral equations over [0, 1] are derived, 
providing an efficient algorithm for the calculation of the functions 
C± in C \ iR* . This algorithm is described in Sec. V. 

To conclude on the cq = 1 case, we note that the functions ( + 
and C- are at the heart of the analytical and numerical calculation 
of the functions X, Y, £x and £y. Obviously, those two auxiliary 
functions are not the only possible steps for the calculation of the 
X- and Y-functions. Different choices can be found in the literature: 
see BusbridgeJ 3 ] Mullikin et al.J 5 ' 7 ! Schultis,! 16 ! Domke,! 17 ! DasJ 18 ! 
and LahozJ 19 ! In these papers, the calculation of the X- and Y- 
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functions amounts to the numerical solution of two Fredholm inte- 
gral equations with regular kernels. The advantage of every proposed 
algorithm lies in the level of simplicity of the investigated integral 
equations, and specially on the level of regularity of their solution. 
Our own scheme, described in Sec. V, is based on the calculation of 
two functions and £_ which are regular-that is defined, continuous 
and differentiable-everywhere on the real axis, including at (con- 
trary to X). Moreover, they are very smooth on the domain [0, 1] of 
the integral equations they satisfy, whatever the values of the albedo 
a G]0, 1[ and optical thickness b g]0, +oo[ are. Finally they satisfy 
the very simple relation (40), which indicates that they are of some 
importance from an analytical point of view. 

A different choice of auxiliary functions was made by Mullikin^ 
forty years ago, leading up to five-digit tables of the X- and Y- 
functionsJ 7 ] Owing to the importance of his work, we have clarified, 
at the end of Sec. V, the link between the functions introduced by 
Mullikin and ours. 

IV. BACK TO THE GENERAL CASE 

The solution to problem (24) for any Co can be deduced from the 
solution C± for cq = 1 using the following relations: 

Co,±(a,M) =Y[M(z)}H(a,z)co(-z) 

+ 2 ^+ ( a ' 6 > z ) ^o,- (a, b,-z)± rjo- (a, b, z)} 

+ 2 C- (a, b, z) [r? ,+ (a, b, —z) =F Vo,+ (a, b, z)] , (57) 

which are valid over C \ iK. The functions ??o,± are defined in the 
whole complex plane by 

r) 0j ±(a,b,z) = ^c T (a,b,+0) 

z r +iac dz' 
+— / H(a,l/z')(C ± na,b,-l/z')c (-l/z')-— -. (58) 

The integral is a principal value at infinity, which can be transformed 
into an integral over [0, 1] using the residue theorem. It is also a 
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Cauchy principal value on the iM-axis, so that the functions ?7o,± are 
sectionally analytic in the complex plane cut along the imaginary 
axis. 

The proof of Eqs. (57) is not given here. It is based on the fact 
that their right-hand sides satisfy Eq. (24), whose unique solution is 

Co,±- 

There is necessarily a link between the functions uo,± and u± 
on the one hand, and between the functions vo t ± and v± on the 
other hand. This is a consequence of the relation (57) connecting the 
functions Co,± and £±. It appears when substituting the expression 
(57) for Co,± into Eqs. (25)-(26), which yields 

uo,±{a,b,z) = \u+(a,b,z)[r]o-(a,b,z)±r]o-(a,b,-z)] 

+ lu-(a,b,z)[r) , + (a,b,z) ^ rf ,+ {a,b,-z)], (59) 

«d i± (a, 6, z) = Y[K(z)]co(-z) ± Y[-^(z)]c (z) exp(tyz) 
+ ±v+(a,b,z)[r)o-(a,b,-z) ± rj ,-(a,b, z)] 
+ ^v-(a,b,z)[f] , + (a,b,-z) Tm,+( a ,b,z)]. (60) 

Replacing the functions uo,± and vo t ± by their above expressions 
in Eqs. (32)-(35), we derive the following solution to Eq. (1): for any 

z G C* 

X (a,b,z) = -[u-(a,b,z)r) 0j+ (a,b,z) + u+ (a,b,z)r) - (a,b,z)}, (61) 

Y (a,b,z) = -[u-(a,b,z)r] 0t+ (a,b-z) - u + (a,b,z)rf -{a,b,-z)], (62) 
and for any z G C \ {— 1, 0} 

co(-z) - %z f 1 X (a, b, v)-^- = Y[X(z)]co(-z) 
* Jo V + z 

+-[v-(a, b, z)r? ,+ (a, b,-z) + v+(a, b, z)rj Q - (a, b, -z)], (63) 
°-z f ' Y (a,b,v)-^- =Y[-X(z)]co{z)exp(b/z) 

i JO V + z 

- 2 l v - ( a > & , z )Vo,+ {a, b, z) - v + (a, b, 2)770 - (a, 6, z)] . (64) 
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Note that Eqs. (60), (63) and (64) are still valid on iR* , provided 
that the Heaviside function Y is extended by 1/2 at 0. 

The above expressions can be confirmed in two particular cases, 
viz. b — ► +oo and Co = 1. When b goes to infinity and the z- variable 
covers the right complex half-plane, Eq. (1) is reduced to an inte- 
gral equation of the form (5), whose solution is given by (II, 50) 
if analyticity is required in the right complex half-plane. This so- 
lution is retrieved by letting b — ► +oo in the expression (61) for 
Xo(a, b, z). We have indeed, using evident notation, (±(a, +oo, z) = 1 
\/z G C, u±(a, +00, z) = H(a,z) for $l(z) > (these functions are 
not defined in the left complex half-plane), v±(a, +oo, z) = 1/H(a, z) 
Vz G C \ { — 1}, and 



rj -(a,+oo,z) 



r) ,+ (a,+oo,z) 
for any z G C. 

Entering these expressions into the general solution (61)-(64) leads, 
as b — ► +00, to 



Xq(q, +oo, z) = H(a, z) 



1 , . z 
2^ + 2ftr 



xco(-lA') 



F (a,+oo,z) = [»(*)>()], 



ff(a,l/0 

[»(z)>0], (66) 
(67) 



co( 



-z) - °-z f 1 X (a, +oo, «)-^- = ¥[»(*)]<*(-*) 
* Jo V + z 



H(a,z) 

%^z I Y (a,+cc,v)- 
^ Jo v 



V + z 

%)+ (o,+oo,-z) (zGC\{-l,0}), (68) 



(if 



+ z 



(z€C\{-l». 



(69) 



Equation (66) coincides with Eq. (II, 50) in the right complex half- 
plane, and Eq. (68) is a generalization of the classical iT-equation 
(II, 31) to any cq(z). The i7-equation is retrieved when Co = 1 since 
7] + {a, +oo, z) = l- Y[-$t{z)]H{a, -z). 
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Another interesting particular case is the one we dealt with in 
Sec. Ill, Co = 1. The functions rj± can easily be calculated from their 
definition (58) with c$ = 1, applying the method of residues to an 
obvious contour in the right complex half-plane. We find the expres- 
sion 

ri ± {a,b,z) = l-Y[-^{z)]H{a,-z)C±(a,b,z), (70) 

which holds on the iR-axis adopting Y(0) = 1/2. Entering this ex- 
pression into Eqs. (61)-(64) with Co = 1, we do find again the solution 
(50)-(53). 

Coming back to the general case, we note that the definition (58) 
of the functions rjo t ± has to be modified for numerical purposes. The 
integral along the imaginary axis can be transformed into an integral 
over [0, 1] using the residue theorem. We postpone to a forthcoming 
paper developing the numerical solution of problem (1) for any free 
term Co, with applications to cq(z) = exp(— rjz) (0 < r < b). The 
main thing is that we know how to calculate the functions ??o,± for 
any given cq(z), as soon as we are acquainted with the functions £± 
over [0, 1]. The latter functions thus yield the solution to Eq. (1) not 
only for the particular free term cq = 1, but for any of them as well. 



V. CALCULATION OF THE C±-FUNCTIONS 

Since the functions C+ and £_ are basic to the resolution of singular 
integral equations of the form (1), it is important to derive from 
their definition (37) an algorithm to compute them. As a matter 
of fact, the integral over the iR-axis in (37) may be transformed 
into an integral over [0, 1] using the residue theorem. We derive two 
uncoupled Fredholm integral equations over the range [0,1], which 
have a unique solution and allow the calculation of the C±-functions 
outside [0, 1]. 

Replace H(a,-l/z') by l/[T(a,l/z')H(a,l/z')] in the integrand 
of Eq. (37) and integrate the function 

/ 1 exp(-bz') , 1_ 

Z T(a, l/z>) H*(a, l/z>) ^ b > i/Z >l + zz> 
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along a classical contour in the left complex half-plane (see details in 
Ref. [20]). One gets the following expressions of (±(a, b, z) in C \ iR: 
- if ze C\{] - l,0[U{-l/fc}} 

C±(o,6,z) = l±M±(o,6): '" 



1 + kz 



• / ( 5 /i7 2 )(a,i;)exp(-6/«)C ± (a,6,i;)^,(71) 
- if « = «€]- 1,0[ 



2/cu 

C±(a,6,u) = l±M±(a,6) 



£ [\ g /H 2 )(a,v)exp(-b/v)C±(a,b,v)^-, (73) 



1 + feu 

T {gT/H 2 ){a, -u) exp(6/n)C±(a, 6, -«) 

± ^}(5/^ 2 )(a,«)exp(-6/t;)C ± (a,6^)^, (72) 

where the symbol / denotes the Cauchy principal value of an integral, 

- at z = — 1 

2k 

C ± (a,6,-l) = l T M ± (a,6) — 
=F 

- at z = 

C±(a,6,0) = l, (74) 

- and at z = — 1/k 

C±(a, b, -1/k) = 1 T p(a, b)(' ± (a, 6, +l/fc) 

± 2M ±{ aMS(a,k ) -kb + l^±^ ] 

±^j\ g /H 2 )(a,v)eM-b/vK±(a,b,v)-^. (75) 
We have introduced the coefficient k = k(a), as the unique root in 
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]0, 1[ of the transcendental equation T(a, 1/z) = 0, i.e., 

r(a,l/fc) = l-~lni±| = (0<fe<l). (76) 

The derived coefficients R(a,k), S(a,k), q(a,b) and M±(a,b) are 
defined by 

s ^ = 1 - p + a%ti-ey (78) 

q ^ b) = l H ^t+m eM ~ kb) ' (79) 

M±(a, b) = q(a, b)(±(a, b, +l/k). (80) 

In the integrals of Eqs. (71)-(75), the notation (g/H 2 )(a,v) de- 
notes the ratio g(a,v)/H 2 (a,v) with 

9(a -" ) = r>, ») + (,»,). «>*»<+»■ (Bi) 

Furthermore, we note that the derivatives at z = +l/k of the func- 
tions z — > #(a, z) and z — ► C±( a > z) appear in the right-hand side 
of Eq. (75). They can be deduced from the i7-equation (II, 31) and 
Eq. (71) respectively. 

Equation (71) consists of two Fredholm integral equations over 
[0,1] for the functions z — > (±(a,b,z). A difficulty arises since the 
coefficients M±(a, b) contain the unknown values of £± at It is 

overcome by putting z = +l/k into Eq. (71) in view of calculating 
M±(a, b) as 

M ± (a,b) = 

-^L |1± ^' (a/ ^ )(tI ,„)exp(-6/„)C ± („,^) I ^ ; ],(82) 
which we substitute into the right-hand side of Eq. (71) with 5ft(z) > 
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0. This latter equation then reads in the right complex half-plane as 
C±(a, b, z) = 1 ± t — —r — 

± | J\g/H 2 )(a, v) exp(-b/v)C±(a, b, v) 

_z_ ± g (a6) 2fc* 1 _ 
w+z l=F^(a,6) l + 

The functions C± are calculated by first solving numerically these 
two Fredholm integral equations over [0, 1]. Once they are known in 
the range [0, 1], their values in the remainder of the complex plane 
are deduced from Eqs. (71)-(75) with M±(a, b) as given by Eq. (82). 
Although the free term and kernel of the integral Eqs. (83) are com- 
plicated, their solution can be reached easily and accurately: we com- 
puted them with ten decimal places in Ref. [20] , where the functions 
£± were introduced for the first time.' 21 ! This is due to the smooth- 
ness of these functions over the [0, l]-range: for b > 0, they have a 
finite derivative everywhere, unlike the functions H, X and £x which 
have an infinite derivative at z = 0. 

Tables of the functions and will be published in conclusion 
to their detailed study on both analytical and numerical grounds. 
New tables of the functions X, Y, £x and £y will also be given. The 
algorithm to be used is an improved version of the one described in 
this section. It is based on the resolution of two Fredholm integral 
equations much simpler than (83), since they reproduce Eqs. (83) 
with q(a, b) = 0. 

Link with Mullikin's algorithm for the calculation of the X- and Y- 
f unctions 

Mullikin's algorithm is based on analytic continuation techniques 
which lead to the solution of four Fredholm integral equations. l 5,7 ! 
Constants are also involved, which must be computed from the so- 
lution to these four equations. Mullikin's equations are found again 
by letting z — > +cxd in Eqs. (71), which yields 

C±(a, b, +oo) = 1 ± 2M±(a, b) 

±^ (g/H 2 )(a,v)exp(-b/v)( ± (a,b,v)dv, (84) 
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then by multiplying these equations by kz/(l + kz) and subtracting 
the resulting equations from Eqs. (71). The new unknown functions 

f(z) ± g(z) = 1 ± C -^£ eM-Wz) (85) 



satisfy 



f(z)±g(z) = l±exp(-b/z) 



, exp(-6/z) rA , , . 1 , kz 



H(a,z) L ^ +v ' ' ' ' ^T^i'l + kz 
aexp(—b/z) z f 1 



2 H M l + kz J i9/H)M{1 - kv) 

fill 

X [f( V )±g( V )]-— (86) 

for K(z) > 0. On the other hand, we have 

( T (a, b, +oc) = -=^[1 - £(a„ ± A))(a, 6)], (87) 
Vl — a 2 



(«c\[-i,oo, (88) 



' = jVT^h/H)(a, v)(l - kv) (v £ [0, 1[). (89) 



N(a,v) k 

Equations (87) follow from Eqs. (52)-(53) by letting z — > +oo; the 
coefficients ao and A) are the zero-order moments of the X- and 
K-functions, namely 



ao(a, b) = / X(a,b,v)dv, 0o(a, b) = / y(a, 6, 
Jo Jo 



(90) 



The other two equations are deduced from Eqs. (II, 25) and (II, 20) 
respectively, with 7 = g 1 ^ 2 . 
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Equations (86) then become 
f{z)±g(z) = l±exp{-b/z) 

T ^("o ±A))(a, fr)zexp(-6/z)iV(a, -z) 

T exp(-6/z)iV(a,-z)^z y (i/N)(a,v) 

fill 

*lf( V )±g(v)] — . (91) 

This is equation (4.22) in Ref. [5], whose constants C and D are: 

C = ~(ao + Po)(a,b), D = ^(a -f3o)(a,b). (92) 

Mullikin's functions / and g are thus connected to our £±-functions 
by Eqs. (85). From the very fact that the functions (± are regular 
over the range [— 1,+1] (which will be proved in a forthcoming ar- 
ticle), we deduce that the functions / and g are also regular in the 
domain [0, 1] of the integral equations they satisfy. This is true in 
particular on the right-hand side of 0. However, these two functions, 
contrary to and £_ , diverge at the ends of the range [—1,0] where 
their evaluation is required for computing the X- and Y-functions 
over [0,1]. This behavior generates a loss of accuracy in the values 
of X(a, b, u) and Y(a, b, u) as u tends to + or 1~. 

VI. CONCLUSION 

We solved the singular integral equations of the form (1) by first 
treating the case cq = 1, and then the general case. For Co = 1, the 
unique solution to (1) analytic in C* can be written in terms of 
two auxiliary functions ( + and £_ which satisfy the relation (40). 
These functions are calculated by solving two uncoupled Fredholm 
integral equations over [0,1]. The calculation of the functions X,Y, 
£x and £y follows immediately from that of the £±-nmctions: see 
Eqs. (50)-(53), where u± and v± are given by Eqs. (25)-(30) without 
the 0-subscripts. Incidently, the classical relations for the functions 
X, Y, £x an d £y were retrieved. 

We have then proved that the solution to Eq. (1) for any free term 
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c can be reduced to the one to the same equation with cq = 1. This 
fact undoubtedly gives weight to the £±-functions. The transition 
from the general case to the Co = 1 case is provided by the relation 
(57), which allows writing the general solution Xq to problem (1) and 
related functions in the form (61)-(64). The functions ryo,± appearing 
in this solution are defined by Eq. (58), which shows that they only 
depend on the ^-functions for a given cq. 

Our presentation is thus focused on the functions £ + and the 
computation of which is treated in Sec. V. We could have attempted 
to write the general solution Xq to Eq. (1) in terms of the classical 
X- and Y-functions. Actually, this is just what we did in a first 
(unpublished) version of this article. We got the formula 



which is valid in C*. It yields again the solution (II, 50) in a half-space 
since X(a, +oo, z) = H(a, z) and Y(a, +oo, z) = for $l(z) > 0. We 
have also Xq = X for cq = 1, since 



from the residue theorem. However, we stress that neither the X- and 
Y-functions, nor the £x- an d £y-functions are proper intermediate 
steps to solve the radiative transfer equation within a slab. This 
important question will be discussed in a forthcoming paper while 
solving the integral form of the transfer equation. Anticipating this 
discussion, we conclude that the expression (93) of the solution to 
problem (1) is more elegant but less judicious than the expression 
(61) based on the £±-functions. 

Finally, we note that there are several possible generalizations to 





24 



RUTILY, CHEVALLIER AND BERGEAT 



the H- function in a finite slab. The classical approach consists in 
introducing the couple {X,Y), which yields (H,0) in the right com- 
plex half-plane when the optical thickness b goes to infinity. Another 
possibility is provided by the functions 

H + (a,b,z)=H(a,z)( + (a,b,z), (96) 

H-(a,b,z)=H(a,z)C-{a,b,z), (97) 
which verify from Eq. (40) 

T(a,z)^[H+(a,b,z)H-(a,b-z) + H_(a,b,z)H + (a,b,-z)] = 1. (98) 

These two functions coincide everywhere with the ii"-function as 
b — > +oo [since (±(a, +oo, z) = 1], and we retrieve the factorization 
relation T(a, z)H(a, z)H(a, —z) = 1 by letting b — > +oo in Eq. (98). 
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APPENDIX: PROOF OF EQ. (40) 

Equation (40) can been proved by first calculating the product 
( + (a,b,z)C,-(a,b,—z) from the definition (37) of the functions £±. 
Observing that 

1 1 z' i z" 

(l + zz')(l- zz") ~ z' + z"^l + zz' + 1 - zz"*' 

the product of the two integrals along the ilR-axis is transformed into 
a sum of two repeated integrals. This sum is calculated by inverting 
the order of integration in one of the repeated integrals (which is 
allowed when z ^ iW), replacing then the Cauchy integrals over iW. 
by their expression as given by Plemelj's formulae. I 10 ] After some 
simplifications, one obtains 

. , . \ > / 7 x -, * f +i °° H(a,-l/z') . , 
C+( a, b, z) C - (a, b,-z) = l + - / _ -±L-LJ. exp(-ta') 

x[(( + ) + +((-) + ](a,b,l/z')dz' 

H / TTT^expf-fe') 

z' 2 dz' 

x[(C+) + x(C-) + ](aAiA , ) Tr ^72, 

which shows that the product ( + (a,b, z)(-(a,b, —z) has the form 
1 + ip(z), where ip is an odd function of z. Changing z into —z, we de- 
duce that C+( a i b, —z)(-(a, b, z) = 1 + z) = 1 — <p(z). Hence the 
relation (40) for z£C \ iM* . 



